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Abstract 



We consider the entire characteristic functions of order 2 and we prove 
some decomposition theorems in a multidimensional case. We show that 
the lack of zeros of the density function is a necessary but not a sufficient 
(as in the one-dimensional case) condition for a characteristic function to 
be decomposable. We also find some simple sufficient conditions. 
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t> 

^ ; 1 Preliminaries 

m 

Entire characteristic functions of order 2 with finite number of zeros were con- 
sidered by Lukacs. In pQ and [2] he presents theorems related to characteristic 
functions of the form 

y. (p(t) = P{t)exp(A(t)), teR (1) 

where P and A are polynomials and A is of order 2. P is a polynomial of an 
even degree and has a form 

p(t) = f[(i-h(i+±=), 
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where £j and £.,■ are zeros of P. Lukacs has solved decomposition problem for 
characteristic function of the form (1) where 



ip(t) = P(t)exp 



2+2 1 



oH 



te 



(2) 



Then the density function corresponding to the characteristic function (|2|) 
has the form 



m 



i 



Q{x) exp 



x 
2^ 



aV2ir 
where the polynomial Q can be written as 



x e 



2n 



Q(x) = J2 {-i) k \ k <y- k H k (^) 



k=0 

Here every H). is the Hermite polynomial of order k and A& G K for fc = 
1, 2, ..., 2n. It is also clear that the polynomial Q must be non-negative for all x € 
M. In this paper we will call Q the polynomial associated with the characteristic 
function p. We have two possibilities for characteristic function of the form @: 
either ip is indecomposable or it admits a decomposition 



<f(t) = ip 1 {t)ip 2 {t), 



x e 



(3) 



where (pi and p>2 are non-trivial characteristic functions. There are again two 
possibilities in the case where p> of the form ([3]) is decomposable. From Plucihska 
(and Lukacs) theorem (see [S]) we have 



(a) ip has 
99 2 (t) = P(t)exp 



normal factor p>i(t) — exp 



and then 



, where u\ + <j\ 
(b) <p has factors of the form @ i.e. 



ipj(t) = Pj(t)exp 



a 2 or 



2 + 2 



°P 



, J = 1,2, 



where a\ + a\ = a 2 , P(t) = P 1 (t)P 2 (t), degPi > and degP 2 > 0. 
The following theorems can be find in [2]. 

Theorem 1. Suppose that the characteristic function ip of the form {J|) admits 
a non-trivial decomposition. Then its associated polynomial Q has no real zeros 
(seef^ th.7.3.1). 

Theorem 2. Let ip be the entire characteristic function of the form {J|) and 
suppose that the polynomial associated with ip has no real zeros. Then ip has a 
normal factor (see (a) above and [2] th.7.3.2). 



From the above theorems we see that the possibility of decomposition of a 
characteristic function (of the considered type) is equivalent to the statement 
that the associated polynomial has no real zeros. The aim of this paper is to 
prove theorems related to the decomposition of multidimensional characteristic 
functions. We will consider multidimensional polynomial-normal distribution of 
the form 



/m(x) = ^^p 2 ,(x)exp(-i(x-b) T A(x-b)) , xeR d , 

where P21 is non-negative polynomial of order 21, b G R d and A is a nondegen- 
erate positive d x d matrix. As we know it from the theory of Fourier transfor- 
mation (see also Lukacs [2]) the characteristic function of a polynomial- normal 
distribution is a product of some polynomial and the characteristic function of 
the normal distribution defined by the same matrix A and the same vector b. 
We will prove the following theorem 

Theorem 3. Let the characteristic function ip of d- dimensional polynomial- 
normal distribution has a non-trivial decomposition tp — <p\(p2, where cpi, </? 2 are 
characteristic functions. Then its associated polynomial Q has no real zeros. 

Proof. Let fii^ and / 2 ; 2 be the densities corresponding to characteristic func- 
tions if i and if2 respectively.Then by [3J, Theorem 2, we have 



/2U(X) 



Vdet Ai ( 1 T 
— P2i 1 (x)exp --(x-bj Ai(x 



(27T) 



2 



Vdct A 2 . , / 1 



(2nY- 



T 



/2i 2 (x) = — p 2;2 (x)exp -o( x -t>2) A 2 (x-b 2 ) , (4) 



2 



where P2h, P2i 2 are non-negative polynomials determinated by the zeros of ipi 
and tp% respectively, Ai and A 2 are non-degenerate, positive defined d x d 
matrixes and bi,b 2 G M d . From (g]), from the equality <p(t) = (/?i(t)</> 2 (t) and 
from Borel theorem for Fourier transformation we have 

/2/(x) = / / 2/2 (x - y)f2h(y)dy = 

R d 

■ / p 2 ; 2 (x-y)p 2;i (y)exp ( -- (x - y - b 2 ) T A 2 (x - y - b 2 ) J x 



Vdet Ai det A : 



' nT A /,, K 1 \ ,7, r ™ ^ TDid 



(2n) d 



xexpl --(y-bx) A 1 (y-b 1 ))dy, xG 

(see Maurin [3]). 

Let us assume that the polynomial Q takes value zero at a point xq G 



p<i 



Then / 2 z(x ) = 0, so 
Vdet A x dct A 2 



/ P2/ 2 ( x o-y)P2/i(y)exp (-- (x -y-b 2 ) T A 2 (x y b 2 ; 



(27T) 

K" 

x exp f-± (y - bi) T Ax (y - b x )) dy =0. 

The above equality holds only in case where the integrand function is equal zero, 
because it is non-negative and continuous. Then we have 

/\ P2i 2 (x -y)p2h(y) =o. 

ySR d 

So 

A Pih (xo-y) = V p 2il (y) = 0. (5) 

ySR d 

Since the polynomials p 2 i 1 and p 2 ; 2 are not equal zero, the sets 

{y G M d : p 2(2 (x - y) = 0} and {y G K d : p 2h (y) = 0} 

are closed and their interiors in M. d are empty. Then their sum has the empty 
interior (the trivial case of Baire theorem) and the condition ([5]) is not fulfilled. 
Hence one of the polynomials p 2 i x or p 2 i 2 is equal zero - contradiction. 

Finally, the polynomial Q associated with the given characteristic function 
has no real zeros. □ 

Now we present an example which shows that Theorem 2 is not true in 
multidimensional case. 

Example 4. Let (Xi,X 2 ) be the 2-dimensional random variable with the den- 
sity 

i ( \ 



f(xi,x 2 ) = — (xix 2 - l) 2 + x\ exp'j--^ 



1 / x I 1 / 2 

= 2~P4(a;i,a;2)exp^ -- \x x + x. 
Then the characteristic function has the form 

<P(ti,t a ) = I [t\t\ + 2ti*2 - 24 - t\ + 3] exp j-i (tf +tj)\. (6) 

By |3], Theorem 2 we know that if ip has non-trivial decomposition ip = (pi(f2 
then tpi and <p 2 are characteristic functions of polynomial-normal distributions. 
Let us first prove that in this case it is a product of characteristic function of 
polynomial-normal distribution and characteristic function of normal distribu- 
tion. It is a consequence of the fact that the polynomial 

P(t x , t 2 ) = t\t\ + 2tit 2 -2t 2 2 -tl + 3 



can not be presented as a product of two polynomials of degree 2. 
In fact let us write the polynomial P as 

W(h) = t\ (4 - l) + 2<i*2 - 2*| + 3, 

where for fixed t% it is a polynomial of degree 2 in t\. Then 

A = 4 {2t\ - 4*2. + 3) > for * 2 e M. 

So the roots have the form 



-t 2 ± V 2 *! -4*|+3 
* r 

'-2 



<m = Vi ■ (7 » 



If the polynomial P is decomposable there are two possibilities 

1. P(*i,* 2 ) = Qi(*i,* 2 )Q 2 (*i,* 2 ), where 

Qj(ti,h) = Cjtiti + djh + ejt 2 + gj, j = 1, 2. 

The polynomial P is equal zero in (*i,* 2 ) iff one of polynomial Qj ( or both of 
them) is equal zero. Let 

Cj*i*2 + djt\ + Cjt 2 + gj = 0. 

Then for a fixed *2 we have 

* 1 = -^i±f. (8) 

Cj* 2 + Oj 

Formula ([8]) is not of the form ([7]). Indeed homography function has only 
one pole, but roots in (|7J are functions with two poles * 2 = ±1. 

2. Let now 

P(*l,*2)=<9l(*l,*2)Q2(*l,*2), 

where 

Qi(*i,* 2 ) = ai* 2 +ci*i* 2 +diti +ei*2 +.gi 

and 

Q 2 (*i,* 2 ) = M 2 +c 2 *i* 2 +d 2 *i +e 2 *2 +.92- 

For fixed *2 and Qi(*i,* 2 ) = we have 

A = Cj* 2 + 2ci(*i*2 + e*i — 4aiei* 2 — 4aigi, 



-Ci* 2 - di ± ^c 2 *| + 2cidi* 2 + d 2 - 4aiei* 2 - iaiffi 

ti = . 

2ai 



WhenQ 2 (*i,* 2 ) =0 then 



Ml + e 2 * 2 + ff 2 
c 2 i 2 + d 2 



(9) 



(10) 



Functions © and (fTU)) have only one pole whereas the function has two 
poles. Then the polynomial P is indecomposable. We can also obtain this 
result writing P in the form 

P(ti, t 2 ) = (oiti + b\t% + c\t\ti + dxti + eit-2 + gi) (02^! + b 2 t 2 + c 2 tit 2 + d 2 t± + e 2 t 2 + g 2 ) 

and comparing the coefficients in successive powers of variables *i,* 2 . 

Hence, if the characteristic function of the form © is decomposable, then it 
is a product of the characteristic function of a polynomial-normal distribution 
and the characteristic function of a normal distribution. Let 

(f(t u t 2 ) =- [44 + 2*1*2 - 2t| - t\ + 3] exp I --(out? + 2012*1*2 + 022*2) \ x 



x exp <^ -- (t l (1 - an) - 2a 12 ht 2 + (1 - a 22 ) * 2 ; ( ■ 

be such a decomposition. We will show that it is not possible because the 
function 

H>\ (*i , *2) = | [*?*2 + 2 *i*2 - 2*1 - t? + 3] exp |- i (ant? + 2012*1*2 + 032*2) } , 

(11) 
where 

Oii,a 2 2>0, ana22-ai2 > 0, 1-On, l-a 22 > 0, (l-an)(l-a 2 2)-a 2 2 > 0, 

could not be a characteristic function of polynomial-normal distribution. 
The density function for characteristic function of the form dill) is 



h{xi,x 2 ) 



1 

2^ 



1 



[44 + 2*1*2 - 2*1 - *f + 3] 



x exp <,—- {ai\t\ + 2a 12 tit 2 + a 22 t\) \ exp {— it\x\ — it 2 x 2 } dt 2 dt\ = 



27r v /aiia 22 - a\ 2 3 



+ - 



1 



2 V 2 



-x^r 



(011O22 - a\ 2 ) 
6a 12 



r l 4 



2ais 



(an 022 - ai 2 ) 2 
4ai2 



-X 3 Y- 



0-11022 — a 12 



a 22 + 2ai 2 + 



2 
-3 



'12 

022 



, (011022 — a 2 2 ) 2 a 22 \/aiia22 ™ Oi 2 V ana 22 — a 
X 2 / 2 1 



xr- 



-1- 



6a 2 2 



2 / 2 1 2 

a 11^22 — a 12 / a H a 22 — a i2 \ a 22 a ll<222 — &12 



F* 



3a 12 



a 22 



2 2 

a 22 ana 22 — a 12 ana 2 2 — &i 2 



2fli 2 a\ 2 1 
2 "^ 

a 22 «22 a 22 



x.x!>^--(x 2 + r 2 ) 



9 1 2 ^(^y)exp(-^ (x 2 + r 2 ) [ 

27T\/aiia22 - af 2 *> I ^ 



where 



X = 



Xl - ^X 2 
1 022 z 



an 



Y = 



X2 
A22 



Let us see that for an = the coefficient of X 4 disappear and the coefficient of 
X 2 is negative. Then for Y = and for X large enough the polynomial Q has 
negative values at points (X, 0) - contradiction. 

Let us consider the case when 012 7^ 0. Now we prove that the polynomial in 
above density function is non-positive. First suppose that a\i > and substitute 



T 



«12 



011022 ^ «12 



-x, 



a\2 > 0. Then 



Q(X,Y)=T 4 



2y2 



—=T a Y- 

\fOV2 

ai2 1 



T 2 Y 
012 



6 x /ai2 



011^22 ^ 12 



«22 



an 



12 



v ai2 022 012 0,110,22 — o 12 



T^ 



022V a 12 
f 



TY- 



+3- 



3a? 2 



022 



2 2 

022 ana22 — a 12 ana22 — a 12 



,«22 ana22 — a 12 
2ai2 a? 2 1 



Y 2 



0,22 



o , , =:Qi(T,Y). 

a 22 a 22 , 



Let us denote 

Thus 
Qi(T,Y)=p 4 (T,T 

f 



P4(»i,a;2) = (xix 2 - 1) + £ 2 



Y 



IOi2 

1 



a 22 ai2 



1 



6a i2 



»22 0\\022 — 12 Oi2 

O 2 , 3 «?2 



y 2 



012 022 ai2 ana 22 — a 12 
6-1/012 2 4^/012 



011O22 ^ o 12 



/ai2 



"22 



rp2 



TY+ 



(1-2-2 



2 2 

a 22 ana22 — a 12 011022 — 12 



2ai2 a\ 2 1 



022 



a 22 a 2 2 



Let us substitute 



Thus 



Y„ 



I~av2 1 T n 



Tn 



Y, 



1 



and 



Pn I -L m -t-n 



/oi2 n 

Yn_ 
Jorl 



,2' 



Hence 

r\ frr v \ ~ ( ' rr m Yn \ ( 0-12 . . «12 1 6ai 2 \ > 

\ V a 12/ \ a 12 a 22 012 011022-012/ 

1 \ f l x2 



2 i 012 I n- 
a22 ai2 ano-22 — o-yi, 



Gy/au 2 4 v /al^ . , 

3 HF= 1 - ) V«^ (» 



ana22 — a 12 v ai2 fl22 

2 1 A a? 2 \ , 3a? 2 



2 1 ^- 1 — 022 — 2ai2 



2 I J- u,22 ^u-TZ I T / 2 19' 

a 22 ana22 - a 12 V a 22/ («na22 - a \2) 

Then the coefficient B n i of the variable n 2 is equal to 

022 012 1 012 



B n 2 = 

a\2 0,22 0,12 0,\\a22 — a 12 

We will show that B n 2 is negative. 

Let us assume first that B n 2 is not negative. Then 

022 Qi2 1 Q12 

^ T 1 2 • 

ai2 &22 &i2 011022 — a 12 



Successive calculation gives 



2 

2 ^ 1 2 1 a 12°22 

a 2 2 i> 0,22 t «i2 



2 
ail022 — «12 



(because 012 > 0), 

a i2 ( a i2 - 022) ^ a na22 (a-22 + a\ 2 - a| 2 ) (12) 

Since 022 + a\ 2 — a| 2 = 022(1 ~ ^22) + O12 > ancl a n«22 > a 2 2 we have 
ana 2 2 (a 2 2 + a 12 — a 22 ) > a i2 ( a 22 + a?2 — ^22) ■ 
Then by (JTSJ) we obtain 

a 12 — a 22 ^ fl 22 + Ol2 — a 22; 

and therefore 

> a 22 . 

The last inequality is false, because 022 > 0. Then B n 2 < and there exists 
such n £ N, that Qi(T n ,Y n ) < 0. It means that fi(%i,n,X2,n) < at (xx jn , X2, n ) 
corresponding to (T n ,Y n ) - contradiction. 

When a\2 < we shall substitute in the above considerations a\2 by — ai 2 . 
We obtain the same contradiction which means that the function /1 is not a 
density function and tp is indecomposable. 

Then the decomposition theorem, which is true in one-dimensional case is 
not true in d-dimensional case, where d > 1. 



In the next part of this paper we will show that if some condition holds for the 
associated polynomial of characteristic function ip of d-dimensional polynomial- 
normal distribution which has no zeros then the characteristic function is de- 
composable. Let us first prove some auxiliary propositions. 

Proposition 5. Let 

Q(x) = Y a - x " 

\a\<2m 

(where a denotes a multiindex, a — (pc\, . . . , ad) € Zl. and \a\ = Oi\ + . . . + ad) 
be the polynomial of d variables, of degree 2m with real coefficients, positive on 
M. d . Let Q satisfies the following condition 

. Q(xi,...,x d ) n , s 

a := mi = — ; r > 0. 13 

xeiM+ J2 l x "l 

\a\<2m 

Then there exists e > such that, if 

W(x)= Y 6 « xQ ' xeM d 

\a\<2m 

and for every a 

|o Q -b a \ <e, 

then the polynomial W has only positive values on M. d . 
Proof. If x = (xi, ..., x d ) E M. d , then 

W(x) = J2 & " x "= J2 (b a -a a )x a + J2 a ^^ ( 14 ) 

|a|<2m |a|<2m |a|<2rn 

>Q(x)- Y K-a a \\x. a \>Q{*)- I Y e l x "l + e 

\a\<2m \\a\<2r. 

From (fT3"j) we have 

Q(x) 



i+ £ l xC 

|a|<2m 



>o, 



SO 



Q(x)>o 1+ £ |x Q | . (15) 

V \a\<2m J 

From flU]) and (JT5J) we have 

W(x)>a(l + ^ |x«|]- £ [l+ ^ |x Q |J =(o- e )(l+ Y l x ° 

V |a|<2m / V |a|<2m / V \a\<2m 



If e < a, then W(x) > 0, and the theorem is proved. 



□ 



Let us see that for every x £ M. d the following equalities and inequalities 
hold 



i + £N 2m <i+ £ |x»| 

3=1 \a\<2m 

(the left hand side is the addend of the right hand side) ; 



(16) 



i+j2 ix a i = i+]T \x 1 \ ai ..\x d \ ai <i+j2 i+Ew ... U+En 

\a\<2m \a\<2m \a\<2m \ 3=1 J \ j=l 



< 



2»» 



: (i+X>i) + E (i+X>i 

3=1 / \a\<2m \ 3=1 



2m 



[l + E&J i + Ei 



3=1 



where -E^m is the number of elements of the set {a G Z^ : |a| < 2m} ; 



2m 



(17) 



1< 2 7 1+E 



1 2m 



i=i 



< 2 T i+ E |a*l fOT i = i,2,...,d 



fc=l 



and 



2m 



i+Ei 



< 



3=1 



i+EN^+E 



\ j=i j=i ^ fe=i 



i+E 



a;/c 



= (i+d) 2m i+£; 



(18) 



3=1 



From (JTSJ), dTTJ) and JT5) we have 



2)/! 



i+EN 2m ^ 1 + E ix Q i<(i + ^ d m ) i+E^-i ^ 



3 = 1 



|a|<2n 



3=1 



<(i+^ ro )(i+d) 2m i+x: 



1 2m 



3 = 1 



Then we can replace the condition in the above proposition to the one of the 
following equivalent conditions 



. , Q(xi,...,x d ) ^ „ 
b := ml > 0, 

1+ E \ x j\ 
3=1 



(19) 



10 



. Q{x u ...,x d ) 

■■= i n M f d 7 ri^ > °- (2°) 



i+E 



X 



3=1 



Proposition 6. Let Q(x) = ^2 °aX a &e a polynomial of degree 2m with real 

\a\<2m 

coefficients, positive on R d . Then U3]) is equivalent to the following condition 

for every 1 < j < d we have o,j t % m > 0, (21) 

where Hj^m is the coefficient of x^ m in Q. 

Proof. If ej is the j-th vector of the standard base in M. d e 2 — (Sji, Sj2> ■ ■■, djd) > 
where 

1 for i = j, 



Sji "10 for i^ j, 

is the Kronecker's symbol, then for every r £ {0, 1, ..., 2m} the sequence rej £ 
Zl is d— index of length \rej\ < 2m. Let us denote 



T3|) =-> (|2Tj) . Let us see that for every 1 < j < d polynomial of one variable 

2m 



' \ x j e j) ~ 2-~t a i' rX J 1 X J ^ 



r=0 



has positive values on R. Then a^o = ao = Q(0) > and 5j 2m > 0. 
Let a J: 2m = 0. Then 



2m- 1 



x ,— >+oc 



_Q(x^)_ = Um " y _a i ,r. 



lim — X ^ J ^ J/ = li m ^ ' _ — -— - 

7 2m 



+ N 

and 

inf Q }^ = 0. 

xeRd i + EM 2m 

i=i 
The last equality contradicts the condition (fTTJ|). so it contradicts (|13[) . Then we 
must have a,j t 2m > 0. 

(|13p <^= (J2TJ) . Let us write the polynomial Q in the form 

d 
Q{x) ^^a h2m x] m + Yl a « xQ ' x. = (xi,...,x d ) eR d , 

J=l a|<2m 

d 

and let us denote Hxl^ = 5Z l^j'l (Z 1 — norm in R d ). From consideration con- 

i=i 
cerning inequality p7|) follows that for every a £ 7A such that \a\ — r the 



11 



following equality holds 



i+E 



2/u 



< 



3=1 



1 



i+E 



3=1 



d \ 2m (l + ||x|| 1 ) 2m - r ' 



Then taking (p~8|) we have 



2 a a X ° 

|a|<2m 



(1+d)' 



< 



E a Q x c 

a|<2m 



i + E 



(i + NIi 



,2m 



< 



E 2 ll ( 1 + rf ) s 

Z ^ /-i i II || \2m—\a\ v 

|a|<2m I 1 + ll X IU 



3 = 1 



and consistently 



E «» x ° 


E »o x " 


T , |a|<2ra 


|a|<2m 


HxHi-H-oc ^ 2m 
1+ E l*il 
3=1 


11x11^00(1 + HXII!) 2 "* 


us denote 




&o:= 


min a,j^ m - 



Then 



E aj,2mtf m b J2 



i=i 



> 



3=1 



d n i o llxl I, — >+oo 

„2m i i V^ ™2m " M1 



> &0 



i+E^ m i + E 



3=1 



3=1 



(Because in M. d l l — norm H-^ is equivalent l 2m — norm ||-|| 2m j so E ^f™ 

3=1 

(ll X ll 2 m) 2m ~> °°: when H X Hl ~> °°)' 

Let us take R > such that for Hxj^ > i? the following inequalities hold 



E a 

\a\<2m 




d 

i + E 

3=1 





< bo 

3' 



/ j Qj,2rn'Ej 



3 = 1 



i+E 



2, 
" 3 u 



3=1 



12 



Let us denote by Ki(0, R) the closed ball in R d with center at x = and radius 
R with respect to the metric given by the norm 1 1 • 1 1 1 . Then for x £ R d \ K\ (0, R) 



Q(x) 



E a J,2mX 
3 = 1 



2)ii 

i 



E a a X ° 

\a\ <2m 



i+E* 2 

We have 



> 



d 
I. 



E a 3,2mX 



2 m 
3 



i+E4 ra ! + E 



i=i 



i+E 



i=i 



j=i 



E a "X C 

|a|<2m 



i+E 



j=i 



> 



&0 



inf 



Q(x) 



i+E 



> min { — , in f 

3 xeifi(o,i£) 



Q(x) 



i=i 



i + E* 2 

3 = 1 



™2m 



>o, 



because the greatest lower bound of the continuous function 

Q(x) 



ft(x) 



x e 



E* 2 



i=i 



on the compact set Ki(Q,R) is the value of the function at some point x £ 
Ki(0, R), so it is a positive number. It means that the condition (|19p holds, so 
(I2T1) is fulfilled. □ 



Finally from Proposition 5 and Proposition 6 we have 
Theorem 7. Let 



Q(x)= ^ a Q x Q 



|a|<2m 



6e i/ie polynomial with real coefficients of degree 2m, positive on M. d and such 
that the condition \21\) is satisfied. Then there exists e > such that if 

W{x) = Y, 5 « x "' x e M ^ 

a <2m 

is the polynomial with real coefficients fulfilling inequalities 

|o Q - &„| < e, 
/or every a, \a\ < 2m, then the polynomial W takes only positive values on R d . 



We show that property (J21J) is invariant with respect to superposition of the 
polynomial Q with nonsingular linear map of M. d . 

Proposition 8. If a positive polynomial Q £ R [xi, ..., x^] o/ degree 2m satisfies 
condition \21\) and A is an affine isomorphism of Mr, then the polynomial Q a := 
Q o A £ R [xi, ..., Xd] is positive and satisfies condition \21)) . 
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Proof. Let us remind that every affine isomorphism in M. d has the form 

A{x) =F(x) + b, xeK d , 

where F is a linear isomorphism and b £ K d . 

Positivity of the polynomial Qa is clear. Let us assume that b = and 
denote by F _1 the operator norm of the map F _1 with respect to I 1 — norm 
ll-Hj in R d . Then HF -1 ^ > 0. Since Q satisfies (gTJ it satisfies O, too. We 
will show that Qa satisfies (f^CT|) . Let us consider x, y £ M. d such that x =F(y). 
We have Hy^ = H^-^x)^ < H^IL Il x lli and then 

Qa(y) Q(x) 0(x) 



(i + l|y|li) 2m (i + HF-Hx)!!,) 2 " 1 " (l + llF-iiuixiy 2 ™" 



There are two possibilities. 

Case 1: H-F -1 ^ > 1. Then 

l+llF^IIJIxll^HF^I^ + llF- 1 ^ Hxll, = ||F- 1 || 1 (l + ||x|| 1 ) 



and 

Q.A(y) 9(x) > c 

(i + l|y|li) 2m l|i ? - 1 ll? m (i + IM| 1 ) 2m -||F-i||; n ' ' 

where c is the least lower bound from condition (j20]l for Q. Hence 



inr ^— > ^— > 0. 

(i + ||y|li) " HF- 1 ! 1 



l 



Case 2: F" 1 , < 1. Then 



i + IIf- 1 !!, iixii, <i 



Q A (y) Q(x) > 

(i + l|y|| 1 ) 2m "(i + l|x|| 1 ) 2m - ' 



. f Qa{y) . . n 

mf 5— > c > 0. 

y6»-(l + lly|li) " 



and 



This means that 



In both cases the polynomial Qa satisfies condition (|20|) . so it satisfies (|2T|) . too. 
Let us assume now that F — id^d which means that A is the translation by 
vector b = (61, ..., bd) ■ Then from Newton's formula for (y, + 6j) m we have that 
for every 1 < j < d the coefficient of yf m in the polynomial Qa is the same as 
the coefficient of xf m in the polynomial Q. It means that Qa satisfies condition 
(|2Tj) iff Q satisfies this condition. Thesis of the proposition follows from the 
fact that every affine isomorphism is a composition of a linear isomorphism and 
some translation. □ 
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Conclusion 9. In the above proposition we can substitute implication by equiv- 
alence: the polynomial Q satisfies condition I121\) iff Qa satisfies A21\) . 



Proof. If A is an isomorphism then A 1 is an isomorphism, too and Q — 
(Qa) a -i ■ □ 

The following theorem is the main result of this paper. 
Theorem 10. Let the function f be the density on M. d of the form 



/ (*) = 3-P2/(x) exp ( -- (x - b) T A (x - 1 

(2tt)' V 2 

where pi\ is a positive polynomial of degree 21 satisfying i21\) . and A is a symmet- 
ric and positive matrix of dimension d x d. Let (p be the characteristic function 
of this distribution. Then there are d- dimensional random variables Y and Z 
- the first with polynomial-normal distribution and the second with the normal 
distribution such that 

ip = ipyifiz, 
where ipy and ipz are characteristic functions of Y and Z respectively. 

Proof. Let X be a d-dimensional random variable with density/. Then the 
random variable X := X — b has the density of the form 



Vdct A , , x ( 1 T 



/x(x) = — — —P2i(* + b) exp --x Ax 

(2tt)^ V 2 

The polynomial p^i is positive (from Theorem 8) and satisfies condition ([2~Tj) . 
We have 

p(t) = <px(t)exp(ib-t), teR d . 

If ipx = l PY L pz 1 where Y and Z are d-dimensional random variables which 
satisfies condition of the thesis. Then for every t £ M. d 

ip(t) = exp (zb- 1) <p Y (t)<p z (t) =Vr(t) [exp(ib • t) ip z (t)] = (p Y {t)<fi^(t), 

where Z = Z + b is d-dimensional random variable of normal distribution. 
Hence we may assume that the density / from the thesis of our theorem has the 
form 

v / detA , N ( 1 T 



/(x) = 3-p 2 /(x)exp --x Ax 

(2tt)2 V 2 

Let L be the matrix which reduce A to a normalized diagonal matrix (to 
the identity matrix). 

Let X = LXJ. Then U = L _1 X and the density function of U is given by 
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fv (u) = /i-ix (u) = /x (iu) |det L| = 

detA(dctL) 2 ( \ T 

5 Vii (Lu) exp - - (Lu) A (Lu) 

(2tt)* V 2 



= 3 p 2 ;(u)exp --u J (L 1 AL)u = 

(2tt) 5 V 2 

= rP2i(u)exp (--u T u 

(2tt) 3 V 2 

From Proposition 8 the polynomial p-fz (u) = P21 (Lu) is positive and satisfies 
condition (|21[) . If <£[/ = (^y^z' wnere Y nas a polynomial- normal distribution 
and Z has a normal distribution then 

<^x (t) = ^lu (t) = ^u (i*t) = (p Y (L*t) <^ z (L*t) = ^ (*) Vlz (*) . * e Rd 

We know that Y := LY has a polynomial-normal distribution and Z := UL has 
a normal distribution. Hence, if the thesis of the theorem holds for the identity 
matrix / then it holds for any symmetric and positive matrix A. 

We must only consider the case when the density from the thesis is given by 



-p 2 z(x)cxp --VxN 

(2^ \ 2 U ) 



/(x) = 



Let tp be the characteristic function of /. Then it can be written in the form 



Y 



(t)= ]T /M^rexpf-i^A 

\a\<2l \ 8=1 / 



t e 



where /3 a eR (see [2], §7.3). Let 



y(> 



( t )= j^ ^{itr^L^eni 

\a\<2l \ i=l 



t € 



where 9 € (0, 1). Then the inverse Fourier transform of <pg is equal to 



Mx) = 



(2*)' 



— 1 d 

e ^^n^(y) 

q|<2/ j=l 



exp 






where ay denotes the Hermite polynomial of order ay (see [2], §7.3). It is clear 
that 



lim / e (x) = 

9 -> 1_ (2tt) 



e ^ n ^ fe) 

|a|<22 j=l 



exp 



1 d 



/(x). 
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Moreover the coefficients of the polynomial 

d 



>*(*)= E CII^(t)' xeKd ' 



P2i 

|a|<2( j=l 

converge to the coefficients of the polynomial 

d 



H<2J j=1 

Thus, if we write P22 in the form 

p 2 z(x) = J^ a a x. a 

\a\<2l 

and _ 

P2i,*(x) - J2 M*)* a > 

|a|<2i 

then for every a e Zl , |a| < 2Z we have 

lim b a (9) = a a . 
0->i- 

Let e > be a number from Theorem 7 chosen for the polynomial Q = p2( • Then 
there exists <5 > such that if f — 9 < 5, then \b a (9) — a a \ < e for \a\ < 21. If we 
take 9 G (I — <5, 1) then we get nonnegative polynomial j>2Z,0 on M. d . Hence fg is 
a density function on R d and the random variable Y has PNDd distribution. 
Let 

/z( x ) = 1 J exp ( -7777— -7^ Y] x i 

Then 



z(t)=c^(- 1 -(l-9^tj 
V »=i 



and 

</> Y (t) ^ z (t) = ^ e (t) <p z (t) = </> (t) , t G M. d . 

This ends the proof of our theorem. □ 

Acknowledgments. We express our thanks to J. Wesolowski and Z. Jclonck 
for the inspiration and useful discussions during writing this paper. 



17 



References 

[1] Lukacs, E., Characteristic functions, Griffin, London 1970. 

[2] Lukacs, E., Developments in characteristic function theory, Oxford Univer- 
sity Press 1983. 

[3] Maj, M., Pasternak- Winiarski, Z., Composition and decomposition of mul- 
tidimensional polynomial-normal distribution, J. Math. Sci. Univ. Tokyo, 
No.14 (2007), 511-530. 

[4] Maurin, K., Analysis, Part II, PWN - Polish Scientific Publishers, 
Warszawa; D. Reidel Publishing Company, Dordrecht, Boston London 
1980. 

[5] Plucihska, A., Composition and decomposition of polynomial - normal dis- 
tributions. Math. Society, Proc. of "Fourth Hungarian Colloquium on Limit 
Theorems of Probability and Statistics" 2001. 



18 



